Exact ground states with deconfined gapless excitations for the 3 leg spin— 1/2 tube 
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We consider a spin-1/2 tube (a three-leg ladder with periodic boundary conditions) with a Hamil- 
tonian given by two projection operators — one on the triangles, and the other on the square 
plaquettes on the side of the tube — that can be written in terms of Heisenberg and four-spin ring 
exchange interactions. We identify 3 phases: (i) for strongly antiferromagnetic exchange on the 
triangles, an exact ground state with a gapped spectrum can be given as an alternation of spin and 
chirality singlet bonds between nearest triangles; (ii) for ferromagnetic exchange on the triangles, we 
recover the phase of the spin-3/2 Heisenberg chain; (iii) between these two phases, a gapless incom- 
mensurate phase exists. We construct an exact ground state with two deconfined domain walls and 
a gapless excitation spectrum at the quantum phase transition point between the incommensurate 
and dimerized phase. 

PACS numbers: 75.10.Jm, 75.10.Kt . 75.30.Kz 
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The projection operator approach to spin models pro- 
vided significant results on the ground state properties 
of quantum magnets. Examples include the Majumdar- 
Ghosh Hamiltonian [l|, a spin-1/2 antiferromagnetic 
Heisenberg chain where the two exact ground state wave 
functions are given by a product of purely nearest- 
neighbor valence bonds (pairs of S=l/2 spins forming 
a singlet) with a gapped excitation spectrurn, in accor- 
dance with the Lieb-Schultz-Mattis theorem^. The ex- 
act "valence bond solid" ground state in the Afflcck- 
Kennedy-Lieb-Tasaki (AKLT) model[3| with gapped ex- 
citations is an explicit realization of Haldane's conjecture 
for S* = 1 Heisenberg chains [J] . Further examples include 
the two-dimensional Shastry-Sutherland model [5| which 
has been realized in SrCu2 (603)2 [y]. In the pyrochlore 
lattice, Yamashita and Ucda have introduced a model 
with a macroscopically degenerate ground state [7|. In 
all these cases the Hamiltonian is a sum of projection op- 
erators [8| and positive semidcfinitc by construction, so 
that any state that has energy is an exact ground state. 

Here we extend this approach to a model of spin-l/2's 
arranged in a 3 leg tube geometry, given by the 
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Hamiltonian (see Fig. [T]) . The tube has L triangles and 
periodic boundary conditions are assumed. The projec- 
tor Pi = (4Si • Si - 3)/12, where Si = J^^^^i ^ij is the 
spin operator on the i"^ triangle, gives 1 if the triangle 
has a total spin of 3/2, and if the spin is 1/2. The 
projection R^ acts on the squares that are on the surface 
of the tube. We denote Sq — J2(i j)£a ^*j ^^ the sum of 
the spin operators belonging to the a square plaquette, 
then R^ = (Sq • Sq)(Sq • S^ — 2)/24 projects onto the 



subspace of states where the total spin of the plaquette a 
is 2, and gives if the spin is or 1 (i.e. if a pair of spins 
on the a square form a valence bond). The expanded 
Hamiltonian using the spin operators reads 



L 3 



'H — y ^ y ^ { jj^Sij • Sij+i -i- JiSij- • Si4_i^^ 

+ J2 (Sij- • Si+ij-^i -f Sij- • Si_|-ij-_i) 

+ •^RE [(Si J- • Si-|-i^j)(Si^j + i • Si-|_i,j-|_i) 

+ (Si,j • Sij+i)(Si-|-ij • Si+ij+i) 



-(s. 



)(S. 



^1+1.3 + 1 )\^i.3 + l 



^l.j+1 



)]} 



(2) 



where the intra-trianglc J±^ = Sifn/Gj- 2/^^/3, the inter- 
triangle Ji = 5A'n/6 and J2 = 5il'[]/12, and the ring 
exchange interaction Jre = K^/i. We set K^ = 1 in 
the following. 

Spin tubes [9l4l6| are interes ting not only since there 
exist experimental realizations [l7[ but also as they are 
the next step after spin ladders towards two dimension 
(2D). Batista and Trugman, have recently studied the 
Hamiltonian that is a sum of the R^ operators over all 
the squares of the square lattice and shown that a class of 
states consisting of nearest neighbor valence bond cover- 
ings, where each square plaquette shares a valence bond, 
have zero energy, i.e. are exact ground states ISl. Il9|. 
The spin tube is identical to the square laticc with pe- 
riod three wrapped in one direction. 

Exact diagonalization (ED) of the Hamiltonian ([T|) up 
to L = 12 triangles (36 sites) showed for K/\ = an un- 
usual result: for an odd number of triangles it revealed 
a energy spin-1/2 doublet at k — momentum, while 
for even number we find three singlet ground states: one 
at /c = and two at A: = tt, as shown in Fig. [2ja). 
The appearance of energy eigenstates means that they 
have zero projections with all of the Ra operators in the 
Hamiltonian. As there is no static covering of valence 
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FIG. 1. The spin tube with a snapshot of valence bond cover- 
ing from the exact spin-chiraUty dimerized ground state. The 
shaded plaquettes are not satisfied in this particular configu- 
ration, and only quantum resonance with other configurations 
will make all the plaquettes satisfied. 



bonds on the spin-tube where each plaquette is satisfied, 
this points to a feature that was not encountered in pro- 
jection Hamiltonians so far. Even more striking is the 
appearance of the energy ground state for the tubes 
with odd number of triangles, as in this case necessarily 
a spin is left unpaired. In the following, we will show that 
this apparent contradiction is resolved due to the quan- 
tum mechanical nature of the valence bonds, and that 
the existence of the energy ground state in the odd size 
system is related with the fact that /\ a = is a quantum 
critical point with gapless deconfined excitations. 

Let us begin by considering weakly coupled triangles. 
The Pj splits the 8 states of a triangle into a. S ^ 3/2 
quadruplet and two degenerate S = 1/2 dublets. The 
latter set of four states constitute the kernel of the P^ that 
can be classified according to the spin a and chirality r 
degrees of freedom, with an orthonormal basis spanned 
by the lafrf), where 



|a„±l/2,) 



cr,2\ 



cr,3\ 



(3) 



The iKf'"*) are states with a valence bond between sites 
j + 1, j + 2 and an unpaired spin a on site j of a triangle: 

observing the periodic boundary condition on j. For con- 
venience, we use \r) and |/) for the r^ = ±1/2 chirality. 
In the ii'A — ^ +00 limit the low energy space is 
spanned by the spin-1/2 states given by Eq. (|3]) and the 
effective spin-chirality Hamiltonian reads 
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(5) 
where ai are the spin~-l/2, and fj are the chirality 



„-th 



(pseudospin-1/2) raising/lowering operators of the 
triangle. This is a p articular case of the effective model 
studied in Refs. llCll4l4l.l2d.l2ll. where it has been shown 
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FIG. 2. (a) Energy spectrum from exact diagonalization of a 
tube with 10 triangles and Ka = as a function of momen- 
tum along the tube. The empty (filled) symbols denote spin- 
singlet (triplet) states classified according to the irreducible 
representations (IR) of D3. (b) Low energy two domain wall 
excitations in the thermodynamic limit compared to ED spec- 
tra of small clusters with symmetry compatible with the vari- 
ational solution. The thick line below the (shaded) continuum 
is the bound state. 



that the spectrum is gapped. Since the energy between 
two neighboring triangles is if either the spins or the 
chiralities form a singlet, the positive semidefinite H' has 
a twofold degenerate exact ground state of alternating 
spin and chirality singlet bonds 22] , shown in Fig. [Sja) . 
Actually these two states breaking translational invari- 
ance are ground states not only of the effective model 
([5]) , but are exact eigenstates of Eq. ^ for any value of 
K^, and are ground states for K/^ > 0, when expressed 
as a linear superpositions of valence bonds 23| . In 
this wave function each triangle contains a valence bond 
- this makes the P^ projections on the triangles happy. 
The unpaired spins of the triangles form valence bonds 
that connect pairs of neighboring triangles (they map to 
the spin-singlet bonds of the effective model), as shown 
in Fig. [1] between triangles 2 and 3, and 4 and 5. The 
plaquettes between these connected triangles all have a 
valence bond, so the corresponding i?^s give 0. However, 
out of the three plaquettes belonging to a chirality singlet 
(e.g. triangle 1,2 or 3,4 in Fig.[T]) two have valence bond 
and the third one is seemingly not satisfied. The problem 
can be resolved if we realize that the three |ff '"') states 
used in the definition of | a, r) in Eq. 1^ are not indepen- 



dent, namely \iy° 



a,l\ 



a.2\ 



/^'3\ 



0. So the chirality 



singlet between the i"^ and i + l"^ triangles can be ex- 
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FIG. 3. (a) Schematic representation of one of the two ground 
states of H with ahernating spin and chirahty valence bonds 
(a translation by a lattice vector gives the other one). The 
small dots inside the circle represent the individual spins of 
a triangle, thick solid lines denote valence bonds. The arcs 
connecting two circles stand for a chirality singlet bond \rl) — 
\lr). In (b) and (c) we show the relevant domain walls in odd 
length tubes: (b) |^J), a triangle with S = 3/2 (hexagon); (c) 
the ellipse containing three spin-1/2 triangles with chiralities 
\rrr) + \lll) and centered at i + 1 denotes \rij,-^). Both |^^) 
and \rj^) have spin-1/2 degrees of freedom. 
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In this expression 
the plaquettes between the legs j — 1 and 3, and j = 2 
and 3 both have valence bonds explicitly in each term, 
thus they are satisfied. Using the linear dependence, this 



is exactly the same as 









that 



makes the third plaquette explicitly satisfied. Thus we 
have shown that even though there is no valence bond 
covering that satisfies all the plaquettes, a linear combi- 
nation of 'imperfect' coverings still constitutes a ground 
state (i.e. for any plaquette we can choose a basis where a 
valence bond is explicitly present on that plaquette|2J]). 



For tubes of odd length with periodic boundary con- 
ditions the system cannot be covered with alternating 
spin and chirality singlets, yet ED shows a ground state 
with zero energy. A study of finite size wave functions 
disclosed that the following two types of domain-wall ex- 
citation of the spin-chirality singlet wave function are rel- 
evant: one is a S" = 3/2 triangle connected with valence 
bonds to the two sides, the other is a S* = 1/2 triangle 
connected with chirality triplets to the sides, as shown in 
Fig. [31 We use the notation |^f ) and jryf ) respectively, 
where i denotes the position of the domain wall. |^f ) and 
|?7f ) have eigenvalue with all the R in the Hamiltonian, 
except those that include spins on the i^^ triangle. Using 
these states as a variational basis, we get the following 
non- vanishing overlaps in Fourier-space: {^'\^\^~(.\'^'^ = 
5(l-aLCOsfc)/6, {il\'}i\il) = 5(l-aLC0s/s)/18, and 
(^^|H|77fc) = — 5-\/3(cos/i: — cl) /18, while the overlaps 
are (77^ \'nl) = (1-aLCOsfc) and {H \il) = 1, where 
ah = 8/2^ vanishes for L — > 00. It turns out that for 
ivTA = the v^lCfc^o) + l'7fc=o) i^ ^n eigenstate with 
energy, thus an exact ground state. The propagation of 
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FIG. 4. Relevant two domain wall configurations in the spin 
singlet sector for L even. In (a) the small ellipse denotes a 
chirality triplet {\rl) -I- \lr)), that breaks up into two domain 
walls, (c), (e) and (f) are generally given as \^im 



,t„i 



mn^ 



C— 5^^|), where the domain walls can be 
arbitrarily separated, (b) and (d) shows overlapping domain 
walls, (b) is actually ICjCi+i ~ CiCl+i)) in (d) the chirality 
configuration is \Urr) — \Trll). Arrows connect states between 
which the Hamiltonian 'H,ki^=o has a nonzero matrix element, 
the position of the arrows corresponds to the position of R 
relevant in the overlap. 



the domain walls in an infinite system is given by the 



Et (fc) 
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gapless dispersion. This is uncommon for a wave function 
consisting of short range valence bonds, and is also possi- 
bly observed in Ref. [25|. From the Hellmann-Feynman 
theorem we get that dEQ^/dKi^ ~ 3/4, i.e. the energy 
of the domain wall becomes negative for K/^ < 0, lower 
than the energy of the dimerized state. This indicates a 
phase transition at ATa = between the dimerized and 
a gapless phase. 

Let us now return to tubes with even number of spins. 
In this case the relevant excitations are pairs of domain 
walls that originate from promoting a chirality singlet 
bond into a chirality triplet, as shown in Fig. ID In large 
systems the domain walls can propagate independently 
if they are far from each other. When the domain walls 
get close they can 'overlap' spatially, these components 
can again be retrieved from results of ED for small sys- 
tems. Once all the relevant states are identified, the over- 
lap matrix and the matrix of the Hamiltonian between 
these states can be calculated analytically for any finite 
length systems, allowing us to take the L -^ 00 limit 
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FIG. 5. Schematic phase diagram of the modeL 



[2 31 . At k = TT, a particular linear combination, featur- 
ing dcconfined domain walls, has energy independent of 
the system size, and is therefore the third exact ground 
state. The low energy part of the two domain wall spec- 
trum is shown in Fig. [2jb) : the two independently prop- 
agating walls form a continuum, while the interaction 
between the walls leads to a bound state. The bound 
state is present for any value of k, whereas the energy 
gap between the bound state and the continuum is oc k'^ 
near k = and tt, indicating a fine tuned interaction be- 
tween the domain walls. We note that a finite size gap 
is closing as 1/L^ at fc = in accordance with the Lieb- 
Schultz-Mattis theorem, and that the domain walls in the 
triplet spin sector are also gapless. The two domain wall 
spectrum differs from similar spectra in the Majumdar- 
Ghosh and AKLT model: in the first one a bound state 
emerges from a gapped two spinon continuum at finite 
momenta[26[, while in the latter the gapp ed bound state 
is well separated from the continuum[27|, even if the in- 
teractions are long ranged [28J 

Next, let us explore what happens for the i^A < 0, be- 
low the point where the gap closes. Using ED, we find a 
succession of level crossings where the ground state alter- 
nates between the Ai and A2 symmetry as we decrease 
/Ca from 0. If we recall that the symmetry of a chi- 
rality valence bond is yl2, the alternation is related to 
introducing a pair of domain walls at each level crossing. 
This phase is gapless and incommensurate, however its 
precise characterization we leave for a future work. At 
if A ~ —0.23 the alternation terminates before all the L 
possible domain walls are added, and the incommensu- 
rate phase ends with a first order phase transition. 



For A'a ^ —0.23 the model is in the universality class 
of the spin-3/2 Heisenberg model. In the A'a -^ —00 
limit, where the spins on the triangles are ferromag- 
netically coupled, the efi'ective Hamiltonian is the S = 
3/2 Heisenberg model with bilinear and biquadratic ex- 
changes, that are 7/12 and 1/18 in units of ATq, respec- 
tively. EDs on short chains indicate that the low en- 
ergy spectrum is adiabatically connected to the one of 
the spin-3/2 Heisenberg model, and is therefore gapless. 
In the absence of the four-spin ring exchange interaction 
density matrix renormalization group calculations found 
a l'^* order phase transition between t he g apped dimer- 
ized and the gapless 5^ = 3/2 phase 20|, l29|, with no 
signature of the incommensurate phase. 

In conclusion, we have demonstrated that even though 
the nearest neighbor valence bond coverings can not sat- 



isfy all the projection operators in the Hamiltonian, a 
linear superposition of such coverings is a good ground 
state. Furthermore, we have constructed an exact ground 
state at a quantum phase transition point that involves 
short range valence bonds and deconfined domain walls 
with a gapless excitation spectrum. Identifying the rel- 
evant excitations and combining with exact diagonaliza- 
tion studies, we have conjectured the phase diagram of 
the model (Fig. [5]), with three phases as a function of the 
interaction strength in the triangles: a gapped dimerized, 
a gapless incommensurate, and the gapless Luttinger liq- 
uid phase of the spin-3/2 Heisenberg chain. 
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SUPPLEMENTARY MATERIAL 



dimerized ground state 



The explicit form of the dimerized wave function of alternating spin and chirality singlets reads as 

L/2 



{(T2, 0-4, •••,0'l} '' = 1 



(l/2-<T, 



"(I. 



2i-l 






)-|^: 



"f^2i — 2^2 ao,i 



2i-l 



'2i 



'■>) 



(7) 



The |i/2,_i V2t 



^2»-l ^2t ) 



i'2i~~i ~ ^2i"') term is a chirality singlet between triangles 2i — 1 and 2i. The spin degrees of 
freedom of triangles 2i and 2i-\-\ form a singlet bond | fi) — | it) : therefore the spin of triangle 2i + 1 is the opposite 
to that of triangle 2i, and the —1 factor takes care of the antisymmetrization. 

I 



the ground state with 2— domain walls 

Most of the matrix elements of the Hamiltonian be- 
tween 2-domain wall states, shown in Fig. 4 in the paper, 
can be derived from the one domain wall overlaps. As a 
reminder, the one domain wall onsite terms (^^ | H |^^) = 
5/18 (?7ml^l'7m) = 5/6 and the propagation of a do- 
main wall is described by (^^| "H |ry^_|_i) = — 5a/3/36. 
These overlaps are valid for the infinitely long systems 
(L — >■ 00), for finite systems we have other terms scaling 
as 2-^. 

For tubes of even length, the nonzero overlaps in the 
thermodynamic limit for m and n are sufficiently sepa- 
rated are: 

(ft £^ -F^ £^\'H Ift fi _ fi ft\ ^ 5/9 

\SrriSn S77iSn| '^ |S77iSn ^ni^n/ -'I "^i 

(vLvi - vivl\y- Hnvi - vLvl) =- 5/3, 



(ft £i _ a £t I -^ 

\SmSn S?nSn| '^ 



:t „t 



ct „t 



{vlvi - vivl\n \TiUt±i - '7i;^Cl±i) = -5^3/36. 

For \S.tvi - iLvl) m~n is even, while for |^,t^Cn " CiCi) 



and Ivl^rii - vLvl) m - n is odd. 

When the two domain walls get close, i.e they over- 
lap spatially, both the diagonal and off-diagonal matrix 
element may differ from the general case, and they read 



„t „t 


'lm'lm+1 


n 


„t „i 

'lm'lm+1 


'7m'7m+l/ 


= 5/6, 


Ut ft 


^fift 


n 


ft ft . 


^fift \ 


= 5/12, 


'im'lrn+2 


fi „t 
Sm 'lm+2 


n 


?m'/m+2 


" Sm'7m+2) 


= 5/4, 


ft ft 
SmSrn+1 


^fift 


H 


ft „t 
Sm'/m+2 


" ^mVm+2/ 


= -5/18, 






(Cn+i 


n 


C™+i) 


= 5/6, 




(ft f_i 


+1 


fi ft 
SmSm+1 


n 


Cn+i) 


= -5\/6/36 



These states are ortho normal (i.e. the overlap matrix is 
identity) in the L — >■ oo limit. The matrix elements given 
above were used to produce the variational results in Fig 
2(b). 

The third exact ground state of 'Hk^=o is then given 
as 



-L/2 L-2 L/2 

*2d» = V2C(7r)+2V3|^eKl))+3 ^ (-l)'^|eeK0)-%/3* E (-1)^IC'7K0)- E (-i)'^iw(^,o) 

i=4.odd i=2,ovon Z=l,odd 

(8) 
This wave function is actually valid for a system of arbitrary length L. Above we use the Fourier transform with the 
following phase convention 

m \ \ 

7?77 (fc, = XI I ^m^m,m+i " '?m'?I,m+/ ) ^XP ( Jfc ( m + - 

p\' (9) 

iv {k, l) = Y^ \Cait,m+i - dClm+;} cxp [ik { m ^ 



C(fc) = ^\c.ra+\) Q^V [ik{m+- 



2 



for ^^ (A:, /) and 7777 (fc, Z) / is odd and 1 < / < L/2 since ^^ (fc, /) = — ^^ (fc, — /) , while for ^77 (fc, /) / is even and runs 
from 2 to L — 2. jti + 1/2 corresponds to the center of mass of the two domain walls, and k is the total wave number. 
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